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Abstract 

We investigate the computational aspects of the basket CDS pricing with counterparty risk 
under a credit contagion model of multinames. This model enables us to capture the system¬ 
atic volatility increases in the market triggered by a particular bankruptcy. The drawback of 
this problem is its analytical complication due to its path-dependent functional, which bears 
a potential failure in its convergence of numerical approximation under standing assumptions. 

In this paper we find sufficient conditions for the desired convergence of the functionals associ¬ 
ated with a class of path-dependent stochastic differential equations. The main ingredient is to 
identify the weak convergence of the approximated solution to the underlying path-dependent 
stochastic differential equation. 
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1 Introduction 

It is well known that the first-passage time of a drifting Brownian motion crossing a deterministic 
level is an inverse Gaussian random variable, as its running maximal process of a drifting Brownian 
motion can be characterized with the reflection principle and Girsanov theorem. This result has 
been applied, among others, in studying the default time of a Arm with a structural framework 
in which a Arm defaults at the first time when its asset value falls below its liability value. The 
structural model is acclaimed to have strong economic foundation and is one of the most popular 
models used in pricing single-name credit derivatives. 

It is natural to ask whether one can also characterize the joint distribution of first-passage 
times of correlated drifting Brownian motions crossing some deterministic levels. When there are 
two correlated Brownian motions, [13 [IE] And the joint distribution of Arst-passage times, which 
is an inAnite sum of modifled Bessel functions of the flrst kind. Little is known for Arst-passage 
times of three or more correlated Brownian motions. This limitation makes it difAcult to study 
the default times of multiple names with a Arst-passage time structural model and is one of the 
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reasons that intensity-based reduced-form models are used in pricing portfolio credit derivatives. 
After all, intensity models give more analytic tractability but do not provide economic reasons why 
firms default (see [3 [26] for some recent results in this direction and the references therein). 

In this paper we discuss the pricing of a basket CDS with counterparty (CDS writer) risk. 
The default times of names in a reference portfolio and that of the counterparty are modeled by 
a first-passage time structural model and correlated Brownian motions. Furthermore, we include 
contagion risk in the model, in which the default of a name in the reference portfolio causes a jump 
increase of volatility of the counterparty, which increases the default probability of the counterparty. 
The aforementioned model, for instance [Bj, incorporates both counterparty risk and contagion risk, 
whereby it can realistically explain the severe difficulties experienced by some seemingly default- 
remote banks underwriting super senior tranche CD Os (collateralized debt obligations) during the 
financial crisis of 2007-08. 

Since the joint distribution of correlated default times in a structural model is unknown, we 
compute the price of a basket CDS by the Euler scheme for SDEs and the Monte Carlo simulation. 
As an effective computational tool, Euler scheme has been widely adopted in the credit risk com¬ 
puting for its simplicity and robustness. Is it possible that a price computed from these “robust” 
algorithms could be actually over-valued in the trillions dollar market? 

As we will illustrate in Example the answer is affirmative in general. Then, is there a set 
of broad conditions, which ensures the computation going to the correct value? In this paper, we 
aim to reveal the reasons for the possible mispricing, and further to provide a rigorous justification 
on the mispricing scenarios. To illustrate the idea, we take a simplified example which indeed 
motivates our general setting in Section 

Example 1 Let Vq be the firm value process of the CDS writer, and {Vi : i = those 

of k basket reference names. We assume that the default time of the firm i is given by n = 
inf {t > 0 : V) (t) < 1}, i = 0, 1,. .., A: Moreover, we assume the volatility of firm i follows di ■ {1 + 
a{t)), where 

a(t) = Ehi Urj < t) (1) 

is the total default number of reference names by time t. The main interest of this this paper is to 
compute a premium rate c per annum in the form of 

^ E[/i(ro,...,rfc)] 

E[/2(ro,...,rfc)] 

for some discontinuous functions /i and f 2 with domain see also the exact formula Q. 

Given that Euler scheme with step size h in the approximation of the underlying firm values 
and further their default times , denoted by and respectively, one can approximate premium 
rate c by replacing Ti by r/* in the formula Q, denoted by c^. Our question is that whether the 
convergence c^ ^ c holds as h ^ 0? 

Regarding the numerical schemes of SDEs, there have been extensive research on the topics of 
both weak and strong convergences, see 13 El da da ED and the references therein for excellent 
expositions. In particular, the book m introduces a systematic and rigorous treatment to the 
numerical approximation of the various types of SDEs. To the best of our knowledge, essentially all 
the above literatures on numerical SDEs study the convergence at fixed grid-points in Markovian 
settings. Back to Example [T| the volatility is path-dependent and the presence of first-passage 
times requires information on the whole path. 

In this regard, we turn our attention to study associated mappings from Skorohod path space 
D to real numbers. For instance, the first passage time r* can be rewritten by Ti = 7r(V)) where 
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vr : D —)• M is a mapping defined by 7r(x) = inf{t > 0 : x{t) < 1}. Similar idea can be applied to 
rewrite the premium rate by c = E[i^i(y)]/E[F 2 ( 1 ^)]; where each Fi is a functional on M^+^-valued 
path space with its argument V = {Vi)i=o,...,k £ To this end, our work can be clearly 

divided into two steps: The convergence —)■ c shall be true by the continuous mapping theorem, 
if 


(Ql) converges to V in distribution (see Section 2.3.3) , denoted by =i> V] 

(Q2) Fi and F 2 are continuous almost surely at V with respect to Skorohod topology. 


Regarding (Ql), the weak limit theorems for the whole path in non-Markovian setting can 
be found in m, and our approach for the weak convergence R is also closely related to 

m- However, |13) establishes the convergence based on the continuity assumption of coefficient 
functions, while the volatility of Example is not continuous as a mapping on a path space, and 
hence their result can not be directly applied here. The main reason for the discontinuity is due to 
the dependence on the number of defaults a, see also the tangency problem of m- Nonetheless, a 
as a function on a path space is almost surely continuous with respect to the probability induced by 
V, see Example |10| for details. As such, we bravely attempt to show the desired weak convergence 
under almost sure continuity assumption of coefficient functions, i.e. 

(H) As a mapping from path space to M, the function a is continuous under Skorohod 

topology almost surely with respect to EE“^. 


In the above, refers to the probability measure on induced by V, see more explanation 

in Section [ 2 . 3. 3[ Although (H) is enough for our purpose to cover our motivated example, it is still 
inappropriate since the unknown solution V shall not be included in the assumption. This leads to 
Assumption]^ which serves the same role as (H) with the help of Assumption]^ 

To this end, it is inevitable to go through the entire procedure and carefully reexamine all 
the necessary steps in the weak convergence. Eirstly, we show the tightness of the discrete Euler 
processes, and deduce the convergence of approximating processes to a limiting process almost 
surely by the Skorohod representation theorem. Secondly, we claim the continuity of the limiting 
process, which plays a crucial role in the proof, see Remark 13 Einally, we complete the proof by 


showing that the limiting process is the weak solution of the underlying SDE. 

Regarding (Q2), provided the completion of (Ql), we shall show the convergence in distribution 
fi{T 0 i • • • > ^ fiiTO, • • • , Tfc) of (]^. Although the form of fi corresponding to the pricing for¬ 

mula @ is complicated, it’s enough to examine the weak convergence on the following two simple 
quantities by the continuous mapping theorem (CMT): 


1. (The convergence in single name risk) One shall verify /(r/’' > t) I{ti > t) for arbitrary 

i and t. Applying CMT, it’s sufficient to show E(ri = t) = 0 for any t, i.e. all underlying 
firms has zero probability to get into default at a particular time. This is guaranteed by 
non-degeneracy Assumption ]^ 

2. (The convergence in counter-party risk) One shall verify I{tq > r/*) /(tq > rQ for each 
i > 1. Again by CMT, it’s enough to show E{to = Tj} = 0, i.e. there shall be zero probability 
that two companies default simultaneously. It boils down two sufficient conditions in turn: if 
either (a) CDS writer is independent to all reference names; or (b) CDS writer is not perfectly 
correlated to all reference names and the volatilities are all piecewise constants, then the above 
convergence holds. 
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To close up the introduction, our contribution is summarized as follows. We establish the 
convergence for the approximation of CDS pricing, see Theorem As illustrated above, the major 
mathematical difficulty compared to the existing literature stems from the following added features 
in our model: (a) the contagion risk due to the dependence of the volatility in total number of 
defaults; (b) single/counter-party default risk in terms of first-passage times. As of any other weak 
limit results on the path space, our result has to be established by many assumptions, for which 
we add careful explanations why they are needed, see Section |2.3[ As a result, the seemingly 
cumbersome assumptions can be reduced to a simple condition in Example If cTj > 0 for all i 
and Vo is not perfectly correlated to any of {Vj, i = 1,... A:}, then —>• c holds. 

The paper is organized as follows. Section 2 presents the problem setup and our main result, 
which is to cover rather general scenarios than Example Note that some notions used in Intro¬ 
duction are also slightly extended in an obvious way. Section 3 includes the technical proof of the 
main result. Section 4 includes further discussions related to this work. 

2 Main Results 

2.1 Problem Setting 

Let r = T -|- 1 for some positive constant T. Denote by D” the space of cadlag functions, i.e. 
functions that are right continuous with left limits defined on [0,T] taking values in is 

abbreviated as D. Let (D, T^ {J't}, P) be a filtered probability space satisfying the usual conditions, 
on which a standard A: -|- 1 dimensional Brownian motion W = (Wo, Wi,..., Wfc)^ is defined. Here, 
is the transpose of x. Suppose that /i(-, •) : x [0, T] i-)- and cr(-, •) : x [0,T] 

^(k+l)xik+l) 

are nonanticipating in the sense that ;u(x, t) = /u(x(- At),t) and cr{x, t) = cr{x{- At),t) 
for all t > 0 and x G (see m for details). We consider A; -|- 1 companies with the firm value 

processes V := {Vo, Vi ,..., Vk)'^ satisfying 

dV (t) = diag {V (t)) (/i {V, t) dt + a {V, t) dW {t )), t > 0, (3) 

where the constant V (0) is the firm value at t = 0, diag {V (t)) is a (A:-|- 1) x (A:-|- 1) diagonal 
matrix with diagonal elements V) (t), i = 0,1,..., A; and g = (^Oj • • • j and a = {o'ij)o<ij<k 
represent asset appreciation and volatility rates, respectively. The product of a and its transpose 
reflects the covariance between the movements in the asset values of firms, thus playing a critical 
role in determining the dependence structure among the firm values. 

Let </> stand for g and a. We impose the following structure to cp = g,cr: Eor any x G 
4>{x, •) can be decomposed into a continuous deterministic process 4>c{-) and a pure jump process 
4>j{x, ■) as follows: 


4>{x,t) = (j)c{t) + cj)j {x,t) = (l)c{t) + ^ jfix), 0<t<f, 

i=l 

where functions gj{ respectively aj) : x [0, T] i—>■ (respectively x and : 

pfc-i-i ^ measurable and nonanticipating in the sense that (/>j(x, t) = (l)j{x{- At),t) and 

M^{x,t) = M^{x{-At),t). The function Jf*( respectively Jf) : i—)■ ]R^+^( respectively x 

is measurable for i = 1, 2,.... In the above, M'^{x, t) and jf (x) represent the number of jumps of cf) 
up to time t and the jump size for the Ath jump, respectively. Moreover, the processes 1 1 —>■ 4‘{V,t) 
and 1 1 —)■ Af'^{V,t) are {J^t}t>o-adapted. 
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Without loss of generality, let Vq be the firm value process of a CDS writer and Vi, i = 1,2,... ,k, 
be the firm value processes of the companies in the reference portfolio. Each company has an 
exponential default barrier Li in the form of 

Li (t) = KiV^*^, t > Q,i = Q,l,... ,k. 


where 7 * and Ki are nonnegative constants. The default time for company i is defined by 


Ti := inf {t > H ■. Vi {t) < Li (t)} AT, f = 0,1,..., k 

the first time the firm value falls below the default barrier. We denote by {r(j) : i = 1,2,... ,k} 
the order statistics of {rj : i = 1,..., k}, i.e. T(j) is the fth default time among k companies in the 
reference portfolio. 

Note that the definition of the default time r* is truncated by T = T + 1, but not by the 
maturity T only for its convenience. The advantage is that, our work is reduced to Cadlag space 
from infinite time interval to a finite time interval, while we can keep the probability of {rj = T} 
as zero under non-degenerate condition. This feature will be used to show the weak convergence 
involved with l{T(^i) < T), while preserve the structure of the swap rate defined in Q. 

In pricing derivative securities with the structural model, it is normally assumed that the drift 
coefficient /r of E in Q is equal to the risk-free interest rate r in a risk-neutral setting. We do 
not insist that /i equal to r in this paper as the firm value is not a traded asset and cannot be 
hedged with the no-arbitrage and martingale representation argument. The firm value process V 
is a measure that may have close relation with traded assets but is mainly used to define default 
events. All results still hold if ^ is replaced by r. For the sake of simplicity, the risk-free interest 
rate is assumed to be a positive constant r. The extension to stochastic interest rate model can be 


done under the framework of this paper (see Remark 23 for details). 


A basket CDS is an insurance product in which the underlying is a portfolio of defaultable 
companies and the writer (seller) of the zth default CDS promises to pay 1 — 5* to the buyer of the 
insurance at the ith default time if that happens before the maturity time T of the contract, 
in return the buyer of the fth CDS agrees to pay the writer a premium fee at rate Cj per annum 
on each of pre-specified dates {0 < < f 2 < • • • < fm = T} as long as the ith default has not 

occurred. (In fact, the triggering time for the basket CDS does not have to be the default time of a 
company, it can be any predefined event). If the writer defaults before the maturity of the contract 
or the fth default time, then the CDS contract terminates and there are no further cash flows. The 
risk neutral swap rate Cj is given by 


Ci = 


E (1 - <5,) 1 (r(,) < T) 1 (tq > r(,) A T)] 


E 


ET=i e ini) > h) 1 (^0 > tj) 


(4) 


where Atj = tj — tj-i, j = 1,... ,m, to = 0, and 1 (•) is the indicator function which equals 1 if 
an event occurs and 0 otherwise. The values of tq and r(j) are dependent on the realized path of 
the process V. The evaluation of c* in (|^ involves the expectations of path-dependent functionals, 
which may naturally be computed with the Monte Carlo and the Euler approximation method. 

Let the time interval [0, T] be partitioned into N equally spaced subintervals with grid points 
= nh, n = 0,..., N, h = T/N, and let V^, valued in be the Euler approximating process 

for V, defined recursively by 


:= + 


diag (V^) [Anh + - W{t))) , n = 0 , 


.N-1. 


5 






where ;= /j{V^,nh) and cr^ := a{V^,nh). In the rest of the paper, for ease of the notation 
complexity without ambiguity, we retain the notation to denote the piecewise constant 

interpolation of sequences {{Vj^, ci^) : n = 0, 1 ,... A^}, i.e. 

(i) = {t) = o-m [nh, (n + l)h). (5) 

Random variables Vj^ are J^jj'-measurable, where := J-^h is the information available at time 
Since IR is a Brownian motion, without changing their distributions, we may generate {V^ : n = 
1 , 2 ,..., A^} by the recursive formula 

Kf+i = V!: + diag {v!:) + diag {v^i) a^VhZ^+u ( 6 ) 

where Zn, n = 1,... ,N, are independent k + 1 dimensional standard normal variables and Zi are 
independent of the filtration for I > n and n = 1,...,A^ — 1. 

Corresponding to the Euler approximating process V^, the annual swap rate has the same 
form as that in 0 except To, Tj and r(j) are replaced by Tq , r/* and respectively. In this paper 
we explore under what conditions we have 


lim c,^ = Cj. 

h-5>0 


(7) 


2.2 The Main Results 

In this part, we present the main result after several assumptions. 

Assumption 2 With some positive constant K for all x G and 0 < ti,t 2 < T, for f 

|</>c (ti) - (/>c (^2)1 < a: |ti - , 


= b,a 




Jfix) 


< K, 

< K for all i. 


Assumption 3 a{x,t) satisfies the uniform nondegeneracy condition, i.e. a{x,t)a{x,t)'^ > XI for 
all x G and t G [0, T] fo r some A > 0. 

We next define vr : D x D i—>■ M as 


vr (x, 1) := inf {t > 0 : x (t) < / (t)} A T, ( 8 ) 

which is the first time of the cadlag function x hitting the barrier 1. Let C” be the collection 
of continuous functions defined on [0,T] taking values in M”. is abbreviated as C. For i = 
0,1,... ,k, define two disjoint subsets of the space C as 

Cj = {x G C : TT (x, Li) < T and inf {t > tt (x, Lj) : x (t) < Li (t)} = vr (x, L,)} , 

and 

= {x G C : TT (x, Li) > T} . 

Assumption 4 The mappings x i—)• /i(x, •) and x i—)• a{x, ■) are continuous under Skorohod topology 
at the set {x G : Xj G Cj U C\, i = 0,1,... ,k}. 
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Assumption 5 For all t > 0, {a (t))oj = (<7 (i))jo = 0; o,.s. for I < i < k. 

Assumption 6 a is piecewise constant almost surely, i.e. for strictly increasing stopping time 
sequence {9 q,9i, ...} such that 9 q = 0 and lim^^oo = T, the process a is in the following form 

OO 

a (t) = '^ail{9i-i <t < 9i), a.s., 

i=l 

where Oi is a nonsingular {k + 1) x (k + 1) matrix, measurable with respect to i = 1, 2 ,... 


We now state the main result of the paper, and the explanations are immediately followed by 


Section 2.3 A short remark on convergence in distribution is also included in Section 2.3.3 


Theorem 7 Let V be the k + 1 dimensional Tt-adapted continuous process of the form (© and 
be Euler approximating process for V given by ([^. If Assumptions ^ and[^ hold, then 
converges to V in distribution. If, in addition, either Assumption or jfi holds, then limc^ = Cj for 

all i = 1,2,... ,k. 


h^O 


2.3 Discussions on Assumptions with its Application to Example 

In this part, we discuss the above assumptions combined with Example to illustrate the main 
result. 


2.3.1 Discussion on Assumption!^ 

Assumption is imposed to guarantee the existence and uniqueness of a solution E to Q . As one 
shall further note from Assumption]^ 0(x, •), (fci') and 0j(x, •) are all bounded by some constant 
K. (In this paper, A is a generic constant whose value may change at each line). We may relax 
the Holder-1/2 continuity of (f>c by a condition |(/>c(ti) — </>c(^ 2 )| < 9{\ti ~ ^ 2 !)) where (7 is a bounded 
function satisfying g'(O) = 0 and g{s) tends to 0 as s tends to 0. 


2.3.2 Skorohod space 

To discuss other assumptions, we shall first state the Skorohod metric for space and related 
notions, which are adopted by [T]. Define a uniform metric on by 

\\x-y\\= sup |x(t)-y(t)| ,Vx,y G (9) 

ie[0,T] 


Let A denote the class of strictly increasing, continuous mappings of 


0,T 


onto itself. Then the 


function space is equipped with the Skorohod topology with the metric 

dix, y) = inf {IIA — III V llx o A — yll} , Vx, y G 

AeA 


( 10 ) 


Since d{x,y) < ||x — y\\ the convergence in Skorohod topology does not imply the convergence 
in uniform topology. However, if the limit is in then they are equivalent. 


Proposition 8 (m page 124]) Elements Xn converge to a limit x in the Skorohod topology 

if and only if there exist functions A„ in A such that lim„x,i(A„t) = x{t) uniformly in t and 
lim„ Ant = t uniformly in t. Moreover, if x is in then Skorohod convergence implies uniform 

convergence. 
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2.3.3 Convergence in distribution 


The notion of the convergence in distribution (see HI) or weak convergence (see [I^ ) may be defined 
in various ways in different literatures. Since it plays an important role in Theoremj^and throughout 
the paper, we give a short remark for its clarification. The random elements of our interests 
and V of Theorem]^ are the maps from a probability space (0,^”,P) to a Skorohod metric space 
(pfc+i^^) equipped with Skorohod metric d defined in (10), and often write it as 




I-)- 


Y,fc+1 


,d) 


or V^, 1/ : n I—)• in short if the context has no ambiguity on T, P and d. The distribution of the 

random element V refers to Pl/“^, which is indeed the probability measure on the Borel u-algebra 

i.e. 

py-i(yl) = F{u : V{u) G A}, VA G S(D*^+^). 

Py“^ is sometimes called as the law of V, or the probability measure induced by V, or the push 
forward measure. Similarly, P(l/^)“^ is the distribution of V^. 

We say, as /i —)• 0, converges to V in distribution with respect to Skorohod metric d, if the 
measure P(l/^)“^ weakly converges to Pl/“^, denoted by P(l/^)“^ => Pl/“^. It sometimes called 
converges to V weakly with respect to Skorohod metric. By Portmanteau Theorem, one can 
equivalently define the convergence in distribution in any of five different ways provided by Page 
26 of [I]. For instance, a common definition adopted in many references is that, is said to be 
convergent to V in distribution with respect to Skorohod metric d, if K[f{V^)] —>■ K[f{V)] for all 
/ G C'fe(D^'''^), the space of all bounded continuous (w.r.t. metric d) functions on 

One may already note that the convergence in distribution relies on the topology of 
Indeed, C' 6 (D*^'''^) under uniform topology induced by || • || of ([^ is a bigger space than 06(0^+^) 
induced by Skorohod metric d{-,-) of (10). This immediately yields by the definition that the 
convergence in distribution with respect to the uniform topology implies convergence in distribution 
with respect to Skorohod topology. In the rest of the paper, unless it is specified, the space 
is equipped with Skorohod metric d by default, and convergence in distribution means by default 
the convergence in distribution with respect to Skorhod topology. 


2.3.4 Discussion on Assnmption 

Recall the hitting time operator tt of (|^. This notion enables us to treat the default time r* 
as a function on a random process, i.e. r* = 7 r(Vi,Li). However, one can not assume 7 r(-,Li) is 
continuous in D in general from the following example. 

Example 9 7 r(x, 0) is not upper semicontinuous at x £ C given by 

x{t) = \t- 1/2|, 

since lim„ 7 r(x„,,0) = T > 1/2 = 7 r(x, 0 ) where Xn = x + 1/n. See Figure^for illustration. 

One can also adapt the above idea to illustrate the potential issue arising from the numerieal 
computation of the credit risk model. Let’s assume that the firm value Vi follows the deterministic 
curve x + 1, i.e. I/i(t) = x(t) + l = l + |t — 1 / 2 |, and a tradable derivative price is given by c = ¥.[ti\, 
where ti = inf{t : Vi{t) < 1} A T. One can easily see that an approximation c^ by usual Euler 
seheme does not guarantee the convergence c^ ^ c^. Obviously, it is not a realistic example due 
to its zero volatility. Then, can we avoid misprieing from the computation by assuming non-zero 
volatility of the firm value? 
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Figure 1: Illustration of Example 


The nondegenerate condition Assumption is very important throughout the paper. It not 
only implies any two of firms Vi and Vj are not perfectly correlated, but also makes the barriers 
regular to the underlying diffusion. As an immediate consequence, we have x i—)• 7 r(x, Li) continuous 
PI^~^-almost surely for all i’s. 

Another use of Assumption is on the total number of defaults a in Example One can 
rewrite this a in terms of tt by 

k 

ait) = Y, Li) < t) := d(E) (II) 

1=1 

for some mapping a defined on Example implies that the functional a above is not 

continuous with respect to Skorohod topology, and hence it violates the sufficient condition for 


the weak convergence given in Condition C5.1 of m- However, Example |10| below can verify the 
almost sure continuity under Assumption!^ which is enough for our purpose. 

Example 10 The number of defaults a defined in Example^ ean be rewritten by ( |11[ ). Example^ 
implies that a is not continuous in general. However, under Assumption^ a is continuous under 
Skorohod topology almost surely with respect to Indeed, this follows from the following two 

facts: 

1. By Blumenthal 0-1 law (see m), 7 r(-, Lj) is continuous with respect to Skorohod metric almost 

surely in for all i = 0,... ,k, i.e. 

IP(Vi G {x G D : 7r(-, Lj) is continuous at x}) = 1. 

2. Moreover, I{Tr{-,Lj) < t) is also continuous with respect to Skorohod metric almost surely in 
¥VY , the probability induced by Vj. 

2.3.5 Discussion on Assnmption 

13] requires the mapping x i—)• a{x, •) is continuous under Skorohod topology 


Condition C5.1 in 
in while our volatility violates this condition as of Example 10 


In contrast. Assumption 

may be regarded as a continuity requirement of coefficient functions on a smaller space {x G : 
Xj G CJ U C 2 ,i = 0,1,..., A:}, Eigure provides some examples to illustrate the concept of 
and 6 * 2 . The paths in the union set of CJ and are regular with respect to the boundary Lj 
in the sense that once the path touches the boundary Lj the path pushes through the boundary. 
Continued from Example due to the fact P(Ej G C\A Cf) = 1 under Assumption a also 
satisfies Assumption]^ We refer Section 2.2 of [23| for more detailed descriptions. 
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Figure 2: Illustration of disjoint subsets: Cj and C^. xi and X 2 belong to C\ and C 2 , respectively. 
X 3 doesn’t belong to C\ or C\. 


2.3.6 Discussion on Assnmption or 

Assumption [5 or[^ together with Assumption is to ensure the indicator function 1 (tq > r(j) AT) 
is continuous Pl/^-almost surely. Suppose the CDS writer is default free before maturity T, then 
the indicator function 1 (tq > T(j) AT) is a constant one a.s. For such a case, Assumption 5 

or[^is not required in Theorem [?| From the financial perspective, Assumption implies that CDs 
writer and the reference names are independent of each other from exogenous factors (IFo and Wi are 
independent). This is realistic in practice since an important criteria adopted by practitioners for 
choosing an appropriate CDS writer is that the CDS writer has little correlation with the reference 
names. On the other hand, if the CDS writer has non-zero correlation with reference names (but 
not perfectly correlated due to Assumption then one can still have almost sure continuity of 
1 (tq > A T) by requiring the piecewise constant form of a, which is indeed the case in reality. 
Note that the volatility is calibrated in practice form time to time, not continuously. 


2.4 Other Operators Related to the CDS Pricing Formula 

For the later use, we also introduce other related notions here. For any re G N, define 5”' 
M” X {1, 2,.. ., re} i-A M as 

(x,j) := the jth smallest value among {xi}^^^. 

Finally, define Tj : i-A M, i = 1, 2, as 

Ti (x) := (1 - 5 ,) 1 ( 5 " ({vr (x„, , i) < t) 

1 (n (xo, To) > s’" (^{tT (Xn, Ln)}n=l ,*) A T^ , 


and 


F2 (x) := ^ e '■T Atjl (^S’" (^{vr (x^, T„)}^^^, i'j > tj'j 1 (tt (xq, Tq) > tj ). 
i=i 
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With the help of mappings tt, and Fj, i = 1,2, the default time r* of company i, the rth 
default time r(j) of the reference portfolio, and the swap rates c* and can be expressed respectively 
as follows: 


Ti = 7r{Vi,Li), i = 0,l,...,k, 

. nFijv)] nFijv’^)] 

^ E[F2(1/)] ' E[F2{Vh)Y 

The default time for company i is illustrated in Figure 



Figure 3: Illustration of the default time for company i, Ti = tt (Vi, Li). 


3 Proof of Theorem 0 

In this section we prove Theorem through a number of lemmas. 


3.1 Preliminary Estimates 

We discuss some properties of here, which plays a crucial role in the subsequent parts. Let 
z^{t) := [ji]h. The Euler scheme can be rewritten in terms of interpolated process in the following 
way: 

rV{t) 

V^(t) = Vh, = Vo+ diag{V^(s))iJ^{s)ds + / diag{V^(s))a^{s)dW(s). 

Jo Jo 

For convenience, we also denote AI^ := and AM^ := AI^ — E [AE^|Fj^] . 

Lemma 11 If Assumption^ holds, then satisfies, for p > 2 and h> 0, 


E 

sup 

V^(t) 

p 


_0<t<f 




(12) 
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m 

p- 



E 

sup 



< K 

'z\t2)-z\ti) 



^=1x1 







m 

P 


r 

E 

sup 

E 



< KE 



-=[^1 





/ \V’^is)\Pds 

J 


Tt 




, Vti < t2, (14) 


where [x] is the largest integer not greater than x. 

Proof: For any t £ [h,T], h > 0 and p > 2, using Burkholder inequality and Assumption]^ 
estimate pth moment of V^{t) as follows: 

E\V^{t)\P <E(|yo| + I / diag{V^{s))p^{s)ds\ + | / diag{V^{s))a’^{s)dW{s)\)P 

Jo Jo 


we 


<KE{\Vo\P + \ / diag{V^{s))p^{s)ds\P + \ / diag{V^{s))a^{s)dW{s)\P) 

Jo Jo 

I'P’d) 

<K + E \diag{V^{s))p^{s)\Pds + KE[ \diag{V^{s))a^{s)a'^{sfdiag{V'^{s))\ds]P/‘^ 

Jo Jo 


rzHt) 


<K + K E|l/'*(s)|Pds. 

Jo 

From above inequality and Gronwall’s inequality, we have 

E\V^{t)\P < Ke^. (15) 

For the case of t < h, (15) still holds by noting V’^{t) = 14^(0) when 0 < t < h. For any 
0 < ti < ^2 < F, we have 


sup 



P 




Am+l)h 

sup 

/ diag{V^{s))p^{s)ds + 

/ diag{V^{s))a^{s)dW{s) 


J Z^(tl) 

1 Z^{tl) 



p{rn-\-l)h 

p 

r{m+l)h 

< K sup 

/ diag{V^{s))p!^{s)ds 

+ K sup 

/ diag{V^{s))a^{s)dW{s) 


Jz^iH) 


Jzh(ti) 


rz^(t2) 

diag(y^{s))p^{s) 

P 


/ 

ds + K sup 

/ diag{V^{s))a^{s)dW{s) 

! z^{t\) 



J z^{ti) 


Since the integral of the second term is a martingale, applying ( |15[ ), Assumption]^ Burkholder’s 
inequality and Holder inequality gives us 



m 

p 

rz'^{t2) 


P 

r r^^it 2 ) 


2 

E sup 

E 

<KE / 

diag{V^{s)) 

ds + KE 

/ 

diagiy^{s))a^{s) 

ds 


-=[^i 

J z^{ti) 



. J Z^{ti) 


- 



r^\t2) 


p 





<KE / 

diag{V^{s)) 

ds 


(16) 



Jz^iti) 






<K{t2 -ti + h). 
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Setting ti = 0 and t 2 = T, we obtain 


E 

sup \V^{t)f 

<K\Vo\P + KE 

sup 

m 

p- 


-0<t<f 


Lo<m<[f]-l 

n=0 

- 


<K. 


A review of the proof shows that the inequality (16) still holds when the expectation is replaced 
by the conditional expectation. Hence we have (14). □ 


Lemma 12 If Assumption^holds, then j{V^) 0 as /i —)■ 0, where j{x) := supQ<^^<'p \x{t)—x{t )| 

for X € 


Proof: Since boundedness of and implies that 


/^v:i\<K\v!:\{h+Vh\Zn+i\), 


we have 


j{V^) = sup |AH^| < Khsup \V^\ + KVhsup |H^| sup |^n+i| 


jh jjh 

where in above inequality the first and second term are denoted by and II^, respectively. 
The convergence of E[l^] to zero as /i —)• 0 can be obtained from 

E[I^] < KhE[sup\Vj^\] < i^/i(E[sup|V;^|2])i/2 < Kh. 


The last inequality in the above is due to (12). On the other hand, 

3 


E[Il'*] < 


KVh (: 


n 4 \ 4 


Esup|H„|3 j (Esup|Z„+i| 


= A:/i 4 EsupIlCIs 


3 

,4 \ 4 


< A:/i 4 EsupIlCIs 


3 

,4 \ 4 


hK sup \Zn+i\ 

n 

^ELl\Zr 

N 


1 

|4\ 4 


= A:/i 4 ^Esupll/'^lt j (E|Zi|^)4 

—)■ 0 as h —)■ 0 

The conclusion follows from that j{V^) converges to 0 in as /i —)• 0. 


□ 


Remark 13 Lemma 12 is the key result that enables us to show that the weak limit process V of 
{V^'\ is continuous. Skorohod representation theorem allows us to treat the limit ^ V in almost 
sure sense, that is 

limd{V’^,V) = 0, a.s. (17) 

h 

However, it does not imply almost sure limit with uniform topology, i.e. 

lim ||H^ - H|| =0, a.s. (18) 

h 
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may not he true. In other words, 0 does not imply 

\imf{V\t)) = f{V{t)), Vt 

h 


even for a bounded continuous function f, which is useful in characterizing properties ofV. Indeed, 
one shall prove continuity of V (i.e. P{F G = 1) in advance to make use of (18) from 

Proposition We note that the proof on the continuity of V is missing in 
references therein. 


and the related 


3.2 Weak Convergence of Approximating Solutions 

This part shows that the approximating processes converge to V in distribution. A general 

approach to this goal is first to prove tightness of for extracting a weak limit from any 

subsequence, then to apply the Skorohod representation theorem for passing the limit almost surely, 
and finally to characterize the limiting process as the solution of the underlying SDE, provided there 
exists a unique weak solution. The main result of this subsection is Theorem 

Proposition 14 If Assumption^holds, then there exists a unique weak solution to SDE ([^. 

Proof: For the sake of simple presentation, we assume no jump for /i. Taking logarithm, it is 
sufficient to consider the following SDE: 

dV{t) = p{V,t) dt + a{V,t) dW{t), t>0. (19) 


By definition, V can be constructed uniquely by the following steps: 

1. Let Vi{t) = Vq + fi{s)ds + crc{s)dW{s) for t > 0. 

2. Let Ti = inf{t > 0 : AA'^(Vi, t) = 1} A T. 

3. Let V 2 {t) = Ei(t) for t <ti, otherwise. 


V2{t) 


Vi{ti)+ / fr{s)ds + (ads) + Jf{Vi))dW{s) 


for t > Ti. 

4. Let T 2 = inf{t > ri : M^{V 2 , t) = 2} A T. 

5. Repeat above steps to construct Vi, r*, i = 1, 2,... until r* = T. 


According to |25l Theorem 1.6.3], Vi has a unique strong solution from 0 to T, hence ri is well 
defined. Since ri is no greater than T, Vi has a unique solution from 0 to ri. |25[ Theorem 1.6.3] 
also implies that E|Ei(ti)P < K, so V 2 has a unique strong solution from ri to T with initial value 
Vi(ri) at time ri. Therefore, V 2 has a unique strong solution from 0 to T. Since the number of 
jumps of a is finite due to Assumption we can proceed the above procedure inductively until 
Tn = T, where n is some finite number. Then Vn is the unique strong solution of (19) from 0 to T. 
Hence, the uniqueness of weak solution is ensured according to m Theorem 9.1.7]. □ 


Lemma 15 


If Assumption @ holds, then {P(E^) ^ : h > 0} is tight. 


14 




Proof: According to [H Theorem 3.8.6] it is sufficient to verify condition a) in [U Theorem 3.7.2] 
and condition b) in [U Theorem 3.8.6] in order to show : h > 0} is tight. To verify 

condition a) in [U Theorem 3.7.2] we need to verify that for every e and rational t > 0, there exists 
a compact set Tg^t C such that 

h 


where T^^^ := {x G : inf^gr^ t | X — y\ < e}. It is worth noting that T^j D T 

by Ko the bound of E\V^{t)\‘^ from (12). For every e and t > 0, by setting 6 = 
Fe,t = {a: G : |x| < d}, we have 


Denote 
(^)4 and 


mmv\t) G ri,) =1 - supP(F'^(t) i ri,) 

^>0 h>0 

>i-supP(F^(t) 

h>0 

=1 -supP(|I/'*(t)| > 6) 
h>0 


> 1 — sup 
h>0 


E\V>^{t)\^ 


= 1 - e. 


To verify condition b) in [H Theorem 3.8.6], we need to find some positive /3 and a family {7^((5) : 
0 < (5 < 1, all h} of nonnegative random variables satisfying 

E[\V^{t + u)- A - V^{t -v)\^ M)< E[j’^{6)\J^t] (20) 


for 0 < t < T, 0 < tt < (5, and 0 < u < <5 A t; in addition, 

lim supE[ 7 ^(s)] = 0 (21) 

5^0 h 


and 


limsupE[|l/'^(<5)-D'^(0)|^Al] =0. (22) 

5^0 h 


We claim that /3 = 4 and a family of nonnegative random variables {supQ<^<'r |D^(s)|'^[(l + 2{h A 
2(1)] : d > 0,h > 0} satisfy (20), (21) and (22). For (20), since either \V^{t + rt) — or 

|14^(t) — V^{t — 1 ^) 1 ^ is zero when h > 26 due to piecewise constant form of we then have 


E[\V^{t + u)- A l\J^t]{\V’^{t) - V^{t - v)\‘^ A 1) 

=E[|l/'^(t + u)- A l|T)](|D'*(t) - V^{t - u)\‘^ A l)l(/i < 2d) 


<E[ 


sup 


{t+6+h) 


E ^*4 


\H^{h < 26) 


j-z'-yr 

<E[ 

J Z^{t) 


l^'*(s)|^ds|T'i]l(/i < 2d) 


<E[ sup |F'^(s)|^[d + 2(/iA2d)]|j'i]. 

0<s<T 


(21) follows from (12) and Cauchy Schwartz inequality that 


limsupE[ sup |l/^(s)|‘^[d + 2(/i A 2d)]] < lim Aid. 

<5^0 h 0<s<T <i-^o 
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(22) is shown as follows 


limsupE[|y^((i)-1/^(0)!^] =lim sup E[|y^(5) - y'*(0)|^] 

<5-S’0/i>o '5 -!'Oo</i<5 

m 

< lim sup E[ sup I 

< lim sup K{6 + h + h) = Q. (by ([T^) 

o</i<5 


So, {E(y'*)-i : /i > 0} is tight. □ 

The next result is needed in the proof of Theorem nn 

Lemma 16 (Rosenthal’s inequality, /^7]/) Let p > 2 and (Xi)jgp} be a sequence of independent 
random variables such that, for any n G N and any i G {1, 2,..., n}, E(Aj) = 0 and EdAil^*) < oo. 
Then we have 

n n n 

E(| < Cp max( E|Xd^ (E(x2))p/2), 

2 = 1 2 = 1 2=1 

where Cp is some constant depending on p. 

We now state the main theorem of this subsection on weak convergence of the Euler scheme. 

Theorem 17 Let be the Euler approximating process for V. If Assumptions^ and[^ hold, 
then converges weakly to V as h ^ 0. 


Proof: Since : /i > 0} is tight, hence for an arbitrary infinite sequence, there exists a 

sub-sequence that has a weak limit. We denote this sub-sequence again by {E^}, and its limit 
by V. Due to the uniqueness of weak solution (see Proposition 14), it suffices to show that V 
is the weak solution of (§. Tightness of {E(E^) ^ : /i > 0} implies V is in Moreover, 

since j{V^) 


0 from Lemma 


12 


P is a continuous process due to [U Theorem 13.4]. Using 
Skorohod representation (see m Theorem 9.17]) we can find and V in the same probability 
space (D, {T7},E) such that their distributions are the same as those of and__P, respectively, 
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and Proposition P" converges 
Since is uniformly bounded 


and P^ converges to P almost surely. According to Remark 
to P a.s. under both Skorohod topology and uniform topology, 
from below by a positive constant due to i^sumption p, i = 0,1,..., fc, is regular with respect 
to Li (see [23l Proposition A.l]). Hence (j){V^, •) —)• (j){V, •) a.s. for (j) = p, a in Skorohod topology 
as /i —)■ 0 by Assumption]^ Denote by D{V,(fi) the set {t G [0,T] : </>(P, •) is discontinuous at t}. 
Since the number of discontinuities of (f{V, t) is bounded almost surely, we can write [0, T] \D(P, 4>) 
as the finite union of disjoint intervals L, i = 1,... ,K. Since 4>{V, ■) is Holder-1/2 continuous at 
each li, 4>{V, •) is uniformly continuous at each R. Therefore as h —)• 0, 


,/(P^ •) - </>(D, •) ll/.^O, a.s., l<i<K. 


(23) 


For t G [0, T] denote by 


M (P, t) := V{t)-V (0) - / diag (P (s)) p (P, s) ds. 

Jo 
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Noting that D{V,cj)) is finitely many by Assumptionand t^(0) = y^(0) = Vq, we have 


lim \M{V'^,t) - M{V,t)\ 

h—>0 


h^O h^O 


K 


Y] / diag{V^{s))^i{V^,s) - diag{V{s))^i{V,s)ds 

Jun[o,t] 


< lim 
h^O 


l[0,t] 


diag{V^{s) — V{s))g{V^, s)ds 


+ lim 
h-s>0 


K 


V/ diag{V{s))[g{V'^,s) - g{V,s)]ds 
Jhnlofi 


K 


<A'lim \\V’^-V ||[o,i] H “/^(^r) ll/in[o,t] j diag{V{s))ds 


h-s>0 


i=l 


hn[o,t] 


= 0 . 


Moreover, noting that and have the same law and {M(y^, t) : h > 0} is uniformly integrable 
due to (12), we have 


E[M{V,t)] =E[M{V,t)] = limE[M(l/'*,t)] = lim E[M(y^, t)], 

h—>0 h—^0 


which shows that 


E[M{V,t)] = limE 


'"71=1 




pnh 

J (n—1 


— lim E 
/i-s>0 


Izhit) 


diag 


{n—l)h 


diag V^_i n{V , s) ds 


For 


the last term above, since g is bounded and 14^, = V^{t), we have 


lim E 

h^O 


[ diag (vh^ g{V^,s) ds 

< K lim hE 


Jzh{t) ^ 

/i-s-0 1 

-* 


= 0 . 
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Denote by Mn^ ;= M^{V^ , nh). Using the Euler recursive formula Q and Assumption]^ we have 
E 


E 


^ ^ Z' 77-/l ^ 

(^n-l) (^n-l - Ky'', •s)) 

^ V{n-l)h ^ ^ ) 


< 


E“ 


n=l 


1 

n—1 


V 


L/iJ 


<Khy^Y.Wti +KhY,^ vti K^’"-ACi"i 


n=l 

[si 

n=l 


n=l 


(24) 


sup 

0<t<T 





III 

V^{t) 

+ KhM 

sup 




0<t<T 

n=l 


<Kh^/'^E 

sup 


+ KhE 

sup 



0<t<T 



0<t<T 



which tends to 0 as /i —)■ 0 from ( |l2[ ). Here we have used the fact that J\f^{V^,-) is uniformly 
bounded given Assumption in the last inequality. Therefore, E [M (U,t)] = 0. Using the same 
procedure one can show E [M (V, t) IT's] = M {V, s) for any 0 < s < t, thus M(V, •) is a martingale. 

Denote the Zth component of the vector process X by X;, and let {Xi,Xg) (t) be the cross¬ 
variation of two real processes Xi and Xq up to time t. For t G [0, T] denote by 

Qiq{V,t) := {Vi,Vq) (t) - j {diag{V {s))a{V,s)a{V,sf diag{V ds. 

Again, due to the uniform topology of the convergence of to U, boundedness of and (23), 
and the uniformly integrability of {Qiq{V^,t) : h > 0} from © , we have 

E|Q;g(F,t)| =E|Qz,(U,t)| = limE|Q^g(U^^)| = limE|Q^,(U^^)|. 

Therefore 

Isl 

E|gz,(U,t)| = hm E ^ (Vi% - 


n=l 
rt 


< lim E 
h-s>0 


-I- lim E 
/i-s-O 


-I- lim E 
/i-s-O 


(^diag {s)j cf{V^, s)a{y^, s)^diag {v^ 

til ^ 

[-1 

'-h^ pnh 

E^zEi^^n -1 / {a{V\s)a{V\sf - 
_ 1 J in—I')}! 


n=l 

rt 




(ui) 

) ds 

\ Uv 

/ Iq 


(25) 


18 














































The second and third limit terms in (25) are bounded by 


> n=l 


lim j +KhY,^ [vt-Kn-AK^'^-K^] 


n=l 


+ lim E 
h^O 


''tnKm 


which is zero from a similar argument for showing (24) goes to zero as /i —?■ 0 and (12). For the 

recursive 


E 


<E 


first limit term in (25), using Assumption and the Euler recursive formula ([^ , we have 

r 

n=l 

i-A 

[rI I' fc+l k+l 

E Kn-1 Y^^^^tlUZnA] 


n=l 


+ E 


+ E 


n=l 

[si 


2=1 


2 = 1 


k+l 


<KhE 


E E(‘"5-i)K(<f,ti),4(z«)f - 1| + E (<'ti)i.(<^ti),4(z„)i(Zn)jl} 

n=l i=l l<i^j<k+l 

[ft] k+l 


sup 

0<s<T 


+ E 


K sup |E'^(s)p/i 2 

0<s<T 


E y^AZn)i 


n=l 2=1 


+ E 


[ft] fc+i 

K sup \V\s+h EE[(^-)?-1] 

0<5<T 


22=1 2=1 


-E 


K sup \V^{s+h 

0<s<T 


E E [(ZnAiZ., 

n=l l<ij^j<k+l 




(26) 


The first term of (26) clearly converges to zero. We only show the second term converges to zero. 


The convergence of the other two terms to zero can be shown similarly. Thanks to 0 Theorem 
2-5.7] X]i=i J2iAi{Zn)i is at the order of (^)^/^(log e > 0. Hence 


lim /i 2 
h-s-O 


[ft] k+l 

‘ EE(z»). 

22=1 2=1 


= 0, a.s. 


Therefore once we can justify the exchange of limit and expectation, then the convergence of the 
expectation to zero is shown. To this end, using Lemma p[6| and Young’s inequality, we show uniform 


integrability of {supo<s<r | h 2 


E [ftl y^fc+l/'7 \ 

n=l I-ji=l K^nji 


, h > 0} as follows 




[ft] k+l 

2 

P -| 

E 

K sup \V'^{s)\+^ 

E 'y^AZn)i 

<E 

K sup \V^{s)A 


0<s<T 

n=li=l 


0<s<T 


[ft] k+l 

+ K¥.[hA EE(z«)ii- 

n=li=l 
[ft] k+l [ft] k+l 

<K + max(E E (E E 


2 \ 2 \ 


n=li=l 

5 


n=li=l 


<K + K max(/i°, /i^) 
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which is finite for all small h, here we have used Lemma 16 in the second last inequality above. 
Hence we show (26) converges to zero as h ^ 0. Therefore 


{Vi, Vq) (t) = I (diag {V (s)) a {V, s) a {V, diag [V {s 


TT \T 


Iq 


ds, 


a.s. 


Thanks to Theorem 7.1], there exists a A: + 1 dimensional Brownian motion B such that 
M (y, t) = f* diag (V (s)) a (V, s) dB (s). Therefore, V is the weak solution of Q. □ 


3.3 Properties of First Passage Times 

Lemma 18 For any n € N and j G {1, 2,... ,n}, the mapping x {x,j) from MF to M is 

continuous. 


Proof: The case of n = 1 is trivial. For n > 2 and x = (xi,..., Xn)'^ G K"', functions 

S"' {x, 1) = min (xi,..., x„) and 5” (x, n) = max (xi,..., Xn) 

are obviously continuous. For j G {2,..., n — 1} (a nonempty set only if n > 3), we may decompose 
X into X = (xi,..., Xn-i)^ and x„, and express 5” {x,j) as 

5" (x,i) = min {max (x,j - 1) ,Xn) ,5”"^ {x,j)} ■ 

Given (•,j) is continuous, and 5” (•, 1) and 5"' {■,n) are continuous, we have 5” (•,j) is con¬ 
tinuous for any j G {1, 2,..., n}. □ 

Lemma 19 Given Assumptions^ and\^ hold, then we have 

P (tq = tj) = 0 and P (r(j) = tj) = 0, 

for any i = 1,2,... ,k and j = 1,2,... ,m. 


Proof: Note that firm values Vi, i = 0,1,..., k, satisfy (©• For any j = 1,2,... ,m, we show 
P (tq = tj) = 0 as follows 


P (to = tj) < P {Vo {tj) = Lo {tj)) 


= P Vb(0)exp 


ho {s) - 


{s) 


ds+ bo {s) dBo {s) = Ko exp {'jotj) 


10 


ho (s) - 


bl{s] 


hi 


K, 


0 


ds+ i ho {s) dBo {s) = In 

Jo Fo fUj 


+ 7otj 


(27) 


where 6o('®) = (1) 0; • • ■ j 0)((t(s)cj(s)^)( 1, 0,..., 0)^. Due to the uniform nondegeneracy Assump¬ 
tion and boundedness of a, there exist some positive constants A and A such that 

0 < A < 6o (f) < A for all 0 < f < T, a.s.. 

Thanks to the property discussed in da Appendix A.5], we know the term in ( |27[ ) is zero. Similarly, 
P {p = tj) = 0 for i = 1, 2,..., A:. Finally, P (r(i) = tj) < (a = ^j) = 0- ^ 


Remark 20 In the proof of Lemma 19. we need to show a one-dimensional continuous semi mar¬ 
tingale never hits a given point within finite time almost surely, more specifically see (27). One 
may wonder if uniform nondegeneracy of a can be relaxed to positive definite to show this result. 
The answer is negative. The counter example for one dimensional case is provided in m- 
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Lemma 21 Suppose Assumptions^ is true. Given AssumptionS^orS^ holds, then we have 

P ('To = T(i) A T) = 0, 

for any i = 1,2,... ,k. 

Proof: Since {tq = r(j) A T} C {tq = n AT}, we have 

k 

IP (t-q = pi) AT) < ^ P (to = Ti A T). 

i=l 

Observe that for any i = 1, 2,..., fe, 


P (to = Tj A T) < P (to = n and t* < T) + P (tq = T and n > T) 

< P (2-d random process {(Vq (t), Vi (t))}^ hits the curve (Lo (t), Li (t)) ,0 < t < T) 
= P (process {(Vq (t) — Lq (t) ,Vi{t) — Li (t))}^ hits (0,0) at some time t in [0, T]). 

(28) 


In the above, we use the fact P (to = T and Ti > T) = 0. Similar to the derivation of (27), the 
event {Vi {t) — Li (t) = 0} is equivalent to 


^ bi (s) dBi (s) - In ~ ~ 1 = 0 


Define the processes Vq and V) to be: 

Vj {t) := bj {s) dBj (s) - In - -fjt + (fij (s) - j ds, j = 0, i. 


The differential form of the SDE for D = {Vo,Vi] is as follows 


dV{t) = fii{t) dt + bi (t) Pi (t) dB'{t), 


where 
hi{t) = 




and B' = (^B'q, B'^ is a standard 2-d Brownian motion, pi is an adapted process satisfying 


bi (t) Pi (t) Pi {ty bi (t) = 


(Ti (t) \ ( (Ti it) 


(Tj (t) J \ aj (t) 


and cTj (t) and aj (t) are fth and jth row vector of a (t). Rewriting (28) in terms of the process V, 
target is to show that if Assumption or is true, then 


our 


V (t) = (0,0) for some 0 < t <t) = 0 


(29) 
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Using Girsanov Theorem, we may assume the drift term is zero: 


dV (t) = hi (t) Pi (t) dB' (t) := a (t) dB' (t). 


We first prove (29) is true under Assumption We have pi (t) is an identity matrix for all t > 0 
almost surely. Since we are investigating the first-passage time before T, without loss of generality, 
we assume 

bo (t) = 1 and bi {t) = 1, Vt > T, a.s. 

Then we have 

Vo, Uo) (oo)= oo, (Vi, Vi^ (oo)= oo and (t) = 0, Vt > 0, a.s. 

Define To{t) '■= inf{s : (Vb)Vb)(s) > t} and Ti{t) := inf{s : {Vi,Vi){s) > t}. By the time change 
results for multidimensional continuous local martingales in m Theorem 5.1.9], we have 

B (s) := (Vo (to (s)) , U (f; (s)) ) ^ , s > 0, 

is a standard 2-d Brownian motion under P with initial position ^Vb (0), Vi (0)^ . Since (Vb, bb)(') 

and (Vi, Vi)0 are both continuous and strictly increasing, the inverse maps Tq”^ and T~^ exist and 
are both strictly increasing. It suffices to show 

^B {t) hits origin at some t in ^0, max (T) , T~^ (T)^ ) ~ 

This is true due to max (^) (^)) < oo,a.s., and nonattainability of the origin by the 

2-d Brownian path shown in m Proposition 3.3.22]. 

We next prove (29) is true under Assumption We have ct is a piecewise constant process 
almost surely. In other words, for strictly increasing stopping time sequence {Oq, 0i, ... } such that 
00 = 0 and lim^^oo dn = T + 1 , the process a is in the form of 


d (t) = ^ djl (0i_i <t <di) 


a.s.. 


i=l 


where di is a nonsingular 2x2 dimensional matrix, measurable with respect to Let p be any 

given point in . Partitioning the event in (29), we have 


= p for some 0 < t <T 

= p for some 0i_i <t<9i < T + 1 


s E 

i=l 

= [ F (v{t) = p for some 0. 

V 


V{9i-i) = x\f {y{0i-i) £ dxy 
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But the integrand in the above integral is seen to be zero from the following derivation. 


F I V{t) = p for some 0i_i <t<9i<T + l 


ViOi-i) = X 


' Si-i 


Gi dW (s) = p — V (Oi-i) for some 6i-i <t<6i<T + l 


< sup] 

c 


C dW (s) = p — X for some 9i-i <t<9i<T + l 


'di-1 


= supP [W{t) — W{9i-i) = C (p — x) for some 9i-i <t<9i<T + l 

c 


V{9i.i) = xj 
V(9i.i) = x^ 

V(9i-i) 


= X 


< sup] 
c 

= 0 . 


= C ^ {p — x) for some 0 < t < T + 1 


Here the supremum is taken among any constant 2x2 nonsingular matrix, since di is measurable 
with respect to and constant from 0j_i to 9i. Also, W denotes a 2-dimensional Brownian 

motion starting from zero. The last inequality is using Markovian property of Brownian motion 
and 9i — 9i-i is no greater than T + 1 almost surely. Hence for each i the integral is zero, so is the 
countable summation. □ 

Assumption or 1^ is used to prove ( |29| ). Although it is intuitive to think the two dimensional 
continuous nondegenerate local martingale (without Assumption or should not hit a given 
point within finite time, it seems difficult to prove this rigorously. Since this unsolved question is 
interesting on its own, we describe it in detail below. 

Remark 22 (an open question) Let W be a two dimensional standard Brownian motion with 
respect to a filtered probability space (fi, P,F := {Tt}t>o). Define a continuous local martingale 
Y by 

Y{t) = F(0) + [ a{s)dW{s), t > 0 
Jo 

with H(0) 0. Assume that a is an adapted 2x2 matrix process satisfying 

A||||P < C’"cr{t)G{t)'^< A||^|p, for all (^,t) G x (0, T) with some A, A > 0, a.s. 

The question we have is that whether the following equality 

P{y(f) = (0, 0), for some t G [0, T]} = 0 
is true? We leave this to our future work. 


3.4 Completion of the Proof 

Now we are ready to complete the proof of Theorem Let and Cf 2 be the sets of the 
continuities of functions Fi and F 2 respectively. Due to |22l Lemma 4], with fixed boundary Lj, 
the mapping n {■, Li) is continuous at each x £ Cl U C^- Notice that Fi and F 2 are composition of 
indicator functions, and 7r(-,Lj), and is continuous by Lemma 


18 


Hence we have 


Cf, :={x G C''+^[0,r] ■.Xi£CluCl,i = 0,...,k, and 5''({7r(x„, f) fiT, 

and 7r(xo,Lo) fi 5^({7r(x„, i) AT} C Cf^ 
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and 


Cf 2 :={x G C''+^[0, T]: Xi^C{iJCi,i = k, and 

S^{{TT{xn,Ln}^^i,i) 7 ^ tj, and f{xq,Lq) / tj} C Cp^- 

To apply the mapping theorem (see [U Theorem 2.7]) we need to show P(w : V{uj) G Cp^) = 1 and 
P(a; : V{uj) G Cp^) = 1. It can be seen by 


P(a; : V{u;) G Cp,) 

>P(a; : V{io) G Cp,) 

=F{ViGCiuCli = 0,...,k, and5'^({7r(K,Ln}Li,i) /T, 
and p{Vo, Lo) / S^MVn, Ln)}L„i) A T) 

=F{Vi ^ Cl U C 2 ,i = 0, ■■■ ,k, and r(j) / T, and tq / r(j) A T) 

=P(r(j) / T, and tq 7 ^ A T) (Because P(V) G C] U i = 0,..., /c) = 1 

due to Assumption and [23l Proposition A.l]) 


= 1 . 


(By Lemma 19 and Lemma \li 


Similarly, one can show P(u; : V{ijj) G Cp^) = 1. With AssumptionsandTheorem 17 says that 
=> P as /i —)• 0. Using the mapping theorem, we conclude that 


(v^'^^FiiV), i = l,2, as/i^O. 


Note that T), i = 1,2, are bounded functions, hence {T)(U^) : h > 0} are families of uniformly 
integrable random variables and 


lim E 
h^O 


F, V 


= E[Fi{V)], i = l,2. 


Therefore, lim/i_^o cf = Q 


□ 


Remark 23 (Extension to stochastic interest rate) Suppose the risk-free interest rate r is a Ft- 
adapted continuous process. Let X be a k + 2 dimensional Ft-adapted continuous process which 
represents Vi, i = 0,1,..., A: and r. Define Xi : x N 1 —)> M and I 2 '■ x [0, T] 1 —)> M as 


Xi(x, i) := exp < — 


L 




Xk+i{u)du 


and 


Z 2 {x,t) := exp \ - I Xk+i{u)du 


For G O and x € C we have that Xn converging to x in Skorohod metric implies Xn converging 
to X in uniform metric which implies fg Xn(u)du converging to fQx{u)du for 0 < t < T. This 
shows the mapping from D to M, defined by x ^ fg x(u}du, is continuous at each x G C. Then, 
the two mappings Ti{-,i) and X 2 (-,f) are continuous at {x G : x* G U C^fi = 0,1,..., A:} 

under Skorohod topology. Since the discounting factors exp{— Vudu} and exp{— fg^ Vudu} can 
be rewritten as ZfiXfi) and Z 2 {X,tj) respectively. Theorem^ still holds for the process X under 
the same setting. 
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4 Conclusions 


We have derived the sufficient conditions for the convergence of the approximation of basket CDS 
with counterparty risk under a credit contagion model of multinames by generalizing the known 
weak limit theorems with discontinuous coefhcents under non-Markovian setting. The method 
developed in this paper may be used to study other problems involving running maximal processes 
of correlated Brownian motions, the joint distribution of which is still unknown for dimension 
greater than two. 
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